Fermion Masses in 5*0(10) Models 
Anjan S. Joshipura0 and Ketan M. Patej^ 



Physical Research Laboratory, Navarangpura, Ahmedahad-380 009, India. 



Abstract 

We examine many 50(10) models for their viability or otherwise in explaining all the fermion 
masses and mixing angles. This study is carried out for both supersymmetric and non- 
supersymmetric models and with minimal (10 + 126) and non-minimal (10 + 126 + 120) Higgs 
content. Extensive numerical fits to fermion masses and mixing are carried out in each case as- 
suming dominance of type-II or type-I seesaw mechanism. Required scale of the B — L breaking 
is identified in each case. In supersymmetric case, several sets of data at the GUT scale with or 
without inclusion of finite supersymmetric corrections are used. All models studied provide quite 
good fits if the type-I seesaw mechanism dominates while many fail if the type-II seesaw domi- 
nates. This can be traced to the absence of the h-r unification at the GUT scale in these models. 
The minimal non-supersymmetric model with type-I seesaw dominance gives excellent fits. In the 
presence of a 45// and an intermediate scale, the model can also account for the gauge coupling 
unification making it potentially interesting model for the complete unification. Structure of the 
Yukawa coupling matrices obtained numerically in this specific case is shown to follow from a very 
simple U{1) symmetry and a Froggatt-Nielsen singlet. 
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I. INTRODUCTION 



Grand unified theories (GUTs) which unify strong and electroweak interactions also pro- 
vide a constrained and unified description of the fermion masses and mixing angles. This is 
particularly true in case of theories based on the SO{W) group [l|. All the known fermions 
plus the right handed (RH) neutrino of a given generation are unified into a single spinorial 
16-dimensional representation of the group. As a consequence, in a renormalizable theo- 
ries of this type only three Yukawa coupling matrices Yiq, Yj2q, ^120 and relative strengths 
between them determine six physical mass matrices Mf with / = u,d,l denoting quarks 
and the charged leptons, f = D corresponding to the Dirac mass matrix for neutrinos and 
f = L,R denoting the corresponding Majorana mass matrices for the left handed and right 
handed neutrinos respectively. The labels 10, 126, 120 correspond to three possible Higgs 
representations contained in the product 16 x 16. The presence of all these Higgs fields is not 
necessary and more economical and allowed possibility is to choose only two of them namely, 
10 and 126. The consistent 5*0(10) breaking needs additional Higgs representations - 210 or 
210 -|- 54 - in case of supersymmetric theories and 45 in case of the non-supersymmetric 
50(10) model. 



The minimal supersymmetric 50(10) model with Higgs fields transforming as 10, 126 -|- 



126, 210 has limited numbers of free parameters and is explored in all its details [3|- 
Detailed analysis of fermion spectrum is also presented in case of the non-minimal super- 
symmetric model containing an additional Higgs field in the 120 representation of 50(10) 
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Similar analysis in case of non-supersymmetric 50(10) model is not done and the 
main purpose of the present paper is to provide such an analysis although we also give a 
comprehensive discussion of various numerical fits in the supersymmetric models. 



The non-supersymmetric GUTs do not have built-in explanation of the gauge hierarchy 
problem but they do share several nice features of the supersymmetric models and avoid 
some of the problems associated with the latter. They allow gauge coupling unification 
and can also provide dark matter candidate in the form of axion. In fact, a minimal 
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non-supersymmetric 50(10) model has been revived in recent studies 16|, [iTI]. Earlier 
discussion of gauge coupling unification in such models is given in 18| . The minimal model 
is more economical then the corresponding supersymmetric case as far the choice of Higgs 
representation is concerned and uses only three sets of Higgs fields namely 10//, 126// and 
45//. This choice is argued to lead to successful 50(10) model on two counts: 
(A) Two or three steps breaking of 50(10) to the standard model (SM) is possible through 
the vacuum expectation value (vev) of 45// and 126//. An intermediate scale ~ 10^^ GeV 
allows the gauge coupling unification. This is shown through a detailed analysis using 



two loop renormalization group equations 
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The presence of intermediate scales is also 



welcome from the point of view of explaining neutrino masses. This is unlike the minimal 
supersymmetric model where an intermediate scale required for neutrino masses spoils the 
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gauge coupling unification. 

(B) The minimal model is argued to be complete to the extent that the required pattern 
of the Higgs vacuum expectation values for the gauge symmetry breaking and gauge cou- 
pling unification can emerge from the minimization of the 1-loop corrected Higgs potential. 

The fermion masses in the minimal supersymmetric or non-supersymmetric models arise 
from the following terms 

16i.(riolOH + l%6T26H)16^ (1) 

The above terms represent a part of the superpotential in the supersymmetric case. In the 
non-supersymmetric version, one would need additional Peccei-Quinn like symmetry and 
the above terms would then represent the most general allowed fermion mass terms in the 
model. Two important features of eq.(IT]) are the following. 

1. If the contribution of Yio to the masses of the third generation dominates then one 
gets the b-T unification 

Yb = Yr (2) 

This can lead to large atmospheric mixing if neutrinos obtain their masses from the 
type-II seesaw mechanism 

2. Non-leading contribution to the second generation masses coming from 126^^ imply a 
relation 

= 3ms (3) 
between the muon and the strange quark masses. 

Both these relations are regarded as successful classical predictions of GUTs like SO (10) 
[H or SU{5) a la Georgi-Jarlskog 20|. They are supposed to hold at the grand unification 



scale Mx- The available experimental information does not quite agree with these generic 
predictions at the quantitative level. Extrapolation of the quark masses at the GUT scale 
in the non-supersymmetric theories do not show b-T unification. In supersymmetric the- 
ories, the evolved values of the Yukawa couplings depend on tan /3 and it is found that 
eq.([2]) holds only at some special values of tan /3. Moreover, the presence of supersymmetry 
breaking around weak scale introduces additional finite tan /3 and sparticle mass-dependent 
corrections which need to be included in extrapolation. Eq.([n]) also gets violated in a large 
parameter space at Mx in these theories. Several existing analysis of fermion masses in 
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supersymmetric models |4i-l6l l9 HllUl5[ | are based on simple and somewhat old extrapolation 
of fermion masses presented in 



which does not include finite threshold corrections. As 
far as the non-supersymmetric theories are concerned there has not been any exhaustive 



confrontation of the extrapolated values |22[ of fermion masses and mixing with the simple 
5*0(10) based models. Motivated by this, we address three main issues in this paper. (1) We 
update the existing analysis of fermion spectrum in various supersymmetric models using 



the recent extrapolation of fermion masses and mixing at the GUT scale 23|, |2J| . (2) We 
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provide new fitting of tlie fermion spectrum in various non-supersymmetric models and (3) 
try to understand tlie fitted structure of fermion mass matrices in terms of simple patterns. 
This provides hint into possible flavour structure of fermion spectrum and possible origin of 
large leptonic mixing angles. 

We begin by giving an overview of the existing analysis of fermion masses and mixing. 
Then we discuss our fits in supersymmetric models. Here we consider both the minimal 
and the non-minimal models. Similar analysis is then carried out in Section IV in the 



non-supersymmetric case using the input from [2^. The analysis in these two sections 
demonstrates the viability or otherwise of various 5*0(10) models from the point of obtain- 
ing correct fermion spectrum. In addition, it leads to very interesting structure for fermion 
mass matrices which can be a starting point to uncover the underlying flavour symmetries. 
We identify such structures in Section V and present a summary of all the new results of 
our analysis in Section VI. 



II. OVERVIEW 

Masses of fermions belonging to 16- dimensional spinorial representation of 5*0(10) arise 
from the renormalizable couplings with Higgs fields belonging to 16 x 16 = 10 + 126 + 120 
representations 

16^(ViolOj^ + VT26l26j^ + Vi2ol20^^)16i. (4) 



in a self explanatory notation. The above couplings represent terms of the superpotential 
in case of supersymmetric models and Yukawa interactions in case of non-supersymmetric 
models with Peccei-Quinn symmetry, see latter. (V120), Vio, Yj2q are complex (anti) symmetric 
matrices in generation space. The Higgs fields 10H-|-126iy + 120/^ are forced to be complex in 
case of the supersymmetric theories. The 10 h and 120 h representations of Higgs fields can 
be real in non-supersymmetric models but we shall take them to be complex in this case also 
to allow for the PQ symmetry. Throughout this paper, we shall restrict to the renormalizable 



models. There exists detailed analysis [25[ of fermion masses in non-renormalizable models 
as well. The IOh, 126//, 120// respectively contain 1, 1 and 2 up-like and equal number of 
down like Higgs doublets. It is assumed that only one linear combination of each remains 
light and acquires vacuum expectation value. This results in the following fermionic mass 
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matrices 




Md = H + F + iG , 

= r{H + sF + itu G ), 
Mi = H-3F+ itiG , 
Md = r{H - 3sF + Ud G ), 
Ml = tlF , 

Mr = r-^^F. (5) 

where {G) H, F are complex (anti)symmetric matrices, r, s, ti, tu, to, tl, 'Tr are dimensionless 
complex parameters of which r,rL,rR can be chosen real without lose of generality. The 
effective neutrino mass matrix for three light neutrinos resulting after the seesaw mechanism 
can be written as 

M, = TlF - mMnF-^Ml = M'J + Mi . (6) 

Eqs. flSllBl) describe the most general mass matrices in any renormalizable 5*0(10) models 
and contain a large number of parameters to be of use. Therefore various special cases are 
considered in the literature and we summarize them below. 



A. Minimal Supersymmetric Model 

This model is characterized by the absence of 120h and hence G in eq.([5]). H can be 
diagonalized with real and positive eigenvalues by rotating the original 16-plets in the gen- 
eration space. Hence all the mass matrices are determined by 19 real parameters if only 
type-II or type-I seesaw dominates. These parameters are determined using 18 observable 
quantities. In spite of the number of observables being less than the parameters, not all 
observables can be fitted with required precision due to non-linear nature of eq.(j5|). Eqs.([5]) 
(with G = 0) are fitted to the observed fermion parameters in various papers |^-l6|. The 
most general minimization is performed by Bertolini et al 0] allowing for arbitrary com- 
bination of both the type-I and II seesaw contributions to neutrino masses. The input 
values for quark and the charged lepton masses used in this analysis is taken from 21 
and correspond to tan /3 = 10. The best fits are obtained in a mixed scenario, type-I gives 
slightly worse and type-II scenario is unable to reproduce all the observables within Icr. 
If type-II seesaw dominates then one needs h-r unification at the GUT scale in order to 
reproduce large atmospheric mixing angle. In contrast, the extrapolated values used in the 
analysis do not show complete h-r unification. This results in a poor fit to the atmospheric 
mixing angle at the minimum. Threshold effects can play important role in achieving the h-r 



unification and improves the fit to fermion masses compared to analysis in [6| as we shall see. 
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B. Non-minimal supersymmetric model 



The other case extensively discussed in the literature corresponds to adding a 120-plet of 
Higgs to the minimal model. Fermion masses in models hi,tliis category have been analyzed 
either assuming type-i jgUnl, Q or type-II 

m 
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13, 15 seesaw dominance. In this 



case, the most general model assuming type-II (type-I) dominance has 29 (31) independent 
parameters after rotating to basis with a real and diagonal diagonal H. One needs to make 
additional assumptions in order to reduce the parameter space. Considerable reduction in 
number of parameters is achieved assuming parity symmetry [l3| or equivalently sponta- 
neous CP violation ll|. This leads to Hermitian Dirac mass matrices. In our notation, this 
corresponds to taking all the parameters in eq.(|5]) to be real, see 0, [l3 for details. Such a 
model has only 17 parameters in case of the type-II dominance, two less than in case of the 
minimal model without G but with arbitrary complex parameters. Number of parameters 
can be further reduced by imposing additional discrete symmetry; Z2 10| or /i-r ^ are 



considered in this context. In spite of the reduction in number of parameters the allowed 



fermionic structure is analytically argued 1^, ll3|] to help in reducing tension in obtaining 



correct CP violating phase or fitting the first generation masses. 



Numerical fits depend on whether type-II or type-I seesaw mechanism is used. Compari- 
son of various models in case of the type-II seesaw dominance is made in 15|] . All the models 



in this category give a very good fit to data with a significantly lower than in case of the 
minimal model. The assumption of the type-I dominance leads to better fits compared to 
the type-II case. Moreover, unlike the type-II dominance, one does not need intermediate 
scale [qMi ll. Il4l ] for reproducing the correct neutrino mass scale. This is a welcome feature 



from the point of view of obtaining the gauge coupling unification. All these works are 



based on the use of quark masses derived in 2l| at tan/3 = 10. We shall re-examine the 



non- minimal model with a different set of input which include the finite threshold corrections. 



C. Non-supersymmetric 50(10) models 

One common feature of all fits with type-II seesaw dominated scenarios is the need for 
an intermediate scale Mj ~ 10^^ — 10^^ GeV. This spoils the gauge coupling unification in 
supersjTiimetric theories. In contrast, an intermediate scale in non-supersymmetric frame- 
work helps in achieving the gauge coupling unification. But unlike the supersymmetric case, 
the non-supersymmetric models do not show the b-T unification and thus type-II models 
in this category do not immediately explain the large atmospheric neutrino mixing angles. 
Viability of this scenario can be checked through detailed numerical fits. Unlike the SUSY 
case, there is no systematic and complete three generation analysis of fermion masses within 
non-supersymmetric models. Various issues involved are summarized in a recent paper 
which contains analytic discussion of the simplified two generation case. 
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The most economical possibility for fermion masses and mixing in non-supersymmetric 
model would be to choose a real IOh or 120// and 126// multiplets of Higgs fields. The latter 
is required for neutrino mass generation but by itself, it cannot generate fermion mixing. 



Thus additional 10// or 120// field is also needed. It is argued [26[ that a 126// and a real 
10// cannot fit even two generation case. Thus one needs a complex 10//. Since both the 
real and the imaginary parts of 10// can independently couple to fermions, this would mean 
additional Yukawa couplings. This can be avoided by assigning a Peccei-Quinn (PQ) charge 
to 10//. Consider the following general definition of the PQ symmetry: 



16f ^ e'^GF] 126// -> e^^*"126// 



IOh e-2'°l0//; 120// ^ e-2*"120// ' ^ ^ 

The most general Yukawa couplings allowed by this symmetry once again reduce to eq.©. 
Thus formally both supersymmetric and non-supersymmetric cases look alike. But there 
is an important difference. The renormalization group running of the Yukawa couplings is 
different in these two cases. Moreover the non-supersymmetric case has intermediate scales. 
Thus input values and consequently the resulting fits would be quite different in these two 
cases. 



III. FERMION MASSES IN SUPERSYMMETRIC THEORIES: NUMERICAL 
ANALYSIS 

In this section, we present the numerical analysis of fermion masses and mixing in different 
supersymmetric cases. We use the data in Table (Jll) and define the following function 

i ^ ' 

where the sum i = 1, .., 14 runs over seven mass ratios and four quark mixing parameters 
(given in Tabled!])), ratio of the solar to atmospheric mass squared differences and the solar 



6{r, ) plus the atmospheric (^'23) mixing angles 27|. For the latter we use the values given in 
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These data are fitted by numerically minimizing the function x^- We assume Am^ 



atm 



to be positive corresponding to the normal neutrino mass hierarchy. We also impose the 3a 
upper bound on ^13 while minimizing the x^- Pi denote the theoretical values of observables 
determined by the input expression, eq. ([5]) and Oj are the experimental values extrapolated 
to the GUT scale, cr, denote the errors in O,. 
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A. Minimal Supersymmetric model 



Our input values of the quark masses and mixing angles at the GUT scale are based 
on the analysis in 23| . This uses more precise values of the b and t quark masses and the 
CKM parameters. More importantly, finite threshold corrections induced by sparticles are 



included in this analysis. Analysis in [23[ proceeds in two steps. First, the quark masses and 
mixing angles are determined by fitting the available low energy data and evolving them to 
the supersymmetry braking scale M5. In the second step, finite sparticle induced corrections 
are included and then evolution is performed up to the GUT scale Mx- These corrections 
are expressed in terms of phenomenological parameters 'jd,b,u,t defined below. We reproduce 
their table of values so obtained as Table ([T]) for convenience of the reader. 





A 


B 


C 


D 


CI 


C2 


tan P 


1.3 


10 


38 


50 


38 


38 


lb 














-0.22 


+0.22 


Id 














-0.21 


+0.21 


It 

















-0.44 


y\Mx) 


6+1 


0.48(2) 


0.49(2) 


0.51(3) 


0.51(2) 


0.51(2) 


y\Mx) 


0.0113+[!;J]r 


0.051(2) 


0.23(1) 


0.37(2) 


0.34(3) 


0.34(3) 


y^Mx) 


0.0114(3) 


0.070(3) 


0.32(2) 


0.51(4) 


0.34(2) 


0.34(2) 


Observables 


GUT scale values with propagated uncertainty 


(mu/mc) 


0.0027(6) 


0.0027(6) 


0.0027(6) 


0.0027(6) 


0.0026(6) 


0.0026(6) 


(nid/ms) 


0.051(7) 


0.051(7) 


0.051(7) 


0.051(7) 


0.051(7) 


0.051(7) 




0.0048(2) 


0.0048(2) 


0.0048(2) 


0.0048(2) 


0.0048(2) 


0.0048(2) 


(nic/mt) 


0.0009t[J:J][]Jo6 


0.0025(2) 


0.0024(2) 


0.0023(2) 


0.0023(2) 


0.0023(2) 


{ms/mb) 


0.014(4) 


0.019(2) 


0.017(2) 


0.016(2) 


0.018(2) 


0.010(2) 




0.059(2) 


0.059(2) 


0.054(2) 


0.050(2) 


0.054(2) 


0.054(2) 


(mb/mr) 




0.73(3) 


0.73(3) 


0.73(4) 


1.00(4) 


1.00(4) 


sin 


0.227(1) 


0.227(1) 


0.227(1) 


0.227(1) 


0.227(1) 


0.227(1) 


sin 6*^3 


0.0289tQ;ooo73 


0.0400(14) 


0.0386(14) 


0.0371(13) 


0.0376(19) 


0.0237(18) 


sin ^13 


U.UUZD_o 00045 


0.0036(7) 


0.0035(7) 


0.0033(7) 


0.0034(7) 


0.0021(5) 


Sckm[°] 


56.31 ± 10.24 


56.31 ± 10.24 


56.31 ± 10.22 


56.31 ± 10.22 


56.31 ± 10.27 


56.31 ± 10.25 



TABLE I. The input values of various observables of quark sector and charged lepton masses 
obtained at GUT-scale Mx for various values of tan (3 and threshold corrections 'jt,b.d assuming an 
effective SUSY scale Ms = 500 GeV (see Q for details). 



Column (A)-(D) show the evolved values of quark mass ratios and mixing angles in the 
absence of threshold corrections for various values of tan/3. One clearly sees the absence 
of the b-T unification at the GUT scale except for the low value of tan/3. This changes 
with the inclusion of threshold corrections. These corrections are parameterized by id,u,b,t 
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which are defined in the following manner. The down quark mass matrix is determine by 
the term QY^d^Hd in the minimal supersymmetric standard model. The corresponding term 
QY^d'^H* involving the second doublet H* is not allowed in the superpotential by SUSY but 

it can be radiatively generated after the SUSY breaking. Since tan /3 = 7— , such terms 

\ d/ 

give significant corrections to the tree level values for large tan (3 and should be included in 
evolving fermion masses and mixing from low energy scale to the Mx- The corrected down 



quark matrix is parameterized in |23|, |29| by 

uf^^i + r"* + v^j^f^^ruVcKM)ydiagUR 

where Ul ^ and Vckm are the (diagonal) down quark mass and the CKM matrix before the 
radiative corrections. The loops involving down squark-gaugino generate the second term 
and the loop with up squark-chargino generate the second term. T^m are diagonal in the ap- 
proximation of taking diagonal squark masses in the basis with diagonal quarks. Assuming 
equality of the first two generation squark masses, the diagonal elements Yd = (7^, 7d, 7b) 
correct the down quark masses and r„ = {■ju^lu^lt) correct the CKM matrix in addition. 
The SUSY threshold corrections are included through these parameters and their best fit 
values corresponding to three classical GUT predictions namely eq.([2]) and eq.([3]) and the 
relation = 1 are determined. Last two columns correspond to different values of 7's 
determined this way. Comparison of column C with C1,C2 shows that threshold corrections 
change significantly the b quark mass as well as ^235 ^13- The neutrino masses and mixing 



that we use are the updated low scale values 28|] but the effects of the evolution to Mqut 



on the ratio of the solar to atmospheric mass scale and on the mixing angles are known to 
be small for the normal hierarchical spectrum that we obtain here. 

We now discuss detailed fits to fermion masses and mixing based on the input values 
in Tabled!]). We assume that either the type-I or the type-II seesaw term in the neutrino 
mass matrix dominates and carry out analysis separately in each of these two cases. We can 
rewrite eqQ as follows. 

'3 + s ~ l-s ~ 
'3(1 - s) ~ 1 + 3s ~ 

(Md-Mi), 





= rnir 1 


Md 




Ml 


rLlTlr 

4 


Mr 


4 



■{Md-Mi). (9) 

We have chosen the basis with a diagonal Mi and introduced Md,i = -^Md^i- Thus 

Ml = Diag.{me/mr,mf,/m^, 1) 

Hence all the quantities in the bracket in the above equation depend on the known ratios of 
charged lepton masses. Md is a complex symmetric matrix with 12 real parameters. Since we 
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are fitting the ratios of different mass eigenvalues and mixing angles, the parameter r remains 
free and it can be fixed by m^. (r^) in the case of type-II (type-I) seesaw dominance is 
determined from the atmospheric mass scale. We have total 14 real parameters (12 in 
and complex s) which are fitted over 14 observables. Four unknown observables in lepton 
sector (^^3 and three CP violating phases) get determined at the minimum. Results of 
numerical analysis carried out separately for the type-II and the type-I dominated seesaw 
mechanisms are shown in Table (!IT|) and Table (lllip respectively. Let us comment on the 
results. 





A 


B 


C 


D 


CI 


C2 


Observables 


Pulls obtained for best fit solution 


{mu/rric) 
(mc/mt) 
(nid/ms) 
{ms/mb) 
(me/m^) 
{m^/mr) 
(nib/mr) 


-0.00668428 
0.56521 
-1.21642 
0.112798 
0.0590249 
0.182548 
0.87282 


0.0276825 
0.157569 

-0.891034 
0.440678 
-0.00627804 
0.103214 
2.20829 


0.0259467 
0.0201093 
-0.27664 
0.163272 

0.3944 
0.821485 

2.79368 


0.120767 
0.0730136 
-1.36265 
0.752408 
0.0396087 
0.0192305 
2.34331 


-0.0212532 
0.130288 
-1.04724 
0.884723 
0.0297987 
0.26316 
0.26656 


0.0356043 
0.320944 
-1.57673 
0.789053 
0.0555931 
0.121145 
0.407798 


( ^<ol \ 

\ Am^ J 
sin of 2 
sin 61^3 
sin 0^3 
sin2 e[2 
sin2 0^3 

^CKM [°] 


0.256292 


0.116314 


-0.14908 


0.230056 


0.0188227 


-0.0140039 


0.0730813 
-0.0311676 

1.33502 
0.00836789 

-1.53367 
-0.345931 


0.0702755 
-0.172792 
-0.0354198 
-0.106439 
-4.97038 
-0.163765 


0.0399788 
-0.471738 

0.494732 
-0.599727 

-4.95673 
-0.600814 


0.105989 
-0.0960437 
0.606606 
-0.27881 
-4.70944 
-0.214459 


0.0779176 
-0.757038 
0.890741 
-0.63356 
-2.56294 
-0.650554 


0.127757 
-0.945821 

1.17758 
-0.510182 
-1.84412 
-0.75885 




6.9367 


30.70 


34.52 


30.68 


10.804 


9.3559 


Observables 


Corresponding Predictions at GUT scale 


sin^ 
rLTTir [GeV] 


0.0226508 
19.9399 
337.171 
147.364 
8.37 X 10-1° 


0.0190847 
18.9784 
346.627 
151.912 

6.0 X 10-1° 


0.0206716 
19.5619 
344.795 
146.886 
6.49 X 10-1° 


0.0196974 
11.92 
350.595 
161.702 
6.94 X 10-1° 


0.0239619 
358.789 
12.4786 
194.023 
7.15 X 10-1° 


0.0209208 
1.78569 
349.711 
168.156 

9.1 X 10-1° 



TABLE II. Best fit solutions for fermion masses and mixing obtained assuming the type-II seesaw 
dominance in the minimal SUSY 50(10) model. Pulls of various observables and predictions 
obtained at the minimum are shown for six different data sets. 



• The best fit in the type-II case is obtained at low tan /3 = 1.3 . This case has b-r 
unification and threshold corrections are not very significant. On the other hand, 



10 





A 


B 


C 


D 


CI 


C2 


Observables 


Pulls obtained for best fit solution 


{mu/mc) 
{rric/mt) 
{md/ms) 
(ms/mh) 

{m^Jmr) 
{mb/nir) 


0.0486938 

1.22599 
-0.229546 
-0.932536 
0.0340323 
0.310305 
-0.486477 


-0.180782 
0.130589 
-0.730641 
-0.886438 
0.442759 
-0.526529 
-0.194215 


0.0653101 
0.246294 
0.223201 

-0.977249 
0.103692 
0.881934 

0.0172182 


0.0053847 
0.146932 

-0.748148 
-1.05766 

-0.476364 
0.938701 
-0.34079 


0.0467579 
0.297256 
-2.2904 
0.735548 
0.0649144 
0.705648 
0.789868 


-0.0119661 
0.273346 
-0.689684 
0.000467775 
-0.0648856 
0.0178824 
-0.734937 


( ^^lol \ 

sin 

sin Oin 
sin Of^ 
sin2 0^2 
sin^ 9^23 

5cKM [°] 


0.122267 


-0.10063 


-0.00563647 


-0.120429 


-0.180164 


0.158557 


0.0432634 
-0.281221 

1.37864 
-0.0528379 
-1.22555 
-0.291137 


0.227948 
-0.0401177 
-0.275689 
-0.0598219 
-1.27077 
0.397159 


0.0186715 
-0.167224 
0.926186 
-0.38133 
-1.43475 
-0.350422 


0.084149 
0.0649082 

0.559003 
-0.172148 
0.0548963 
-0.755859 


0.130301 
-0.273222 

1.48675 
-0.746107 

-1.99485 
-0.956628 


0.0922391 
-1.17651 

0.248759 
0.0694831 
-0.946001 

-0.3197 




6.3479 


3.7962 


5.0715 


3.8665 


14.789 


3.4746 


Observables 


Corresponding Predictions at GUT scale 


sin2 
OMNSi J 

a2[°] 


0.0223307 
2.41793 
347.106 
163.759 


0.0194886 
4.52493 
8.42838 
191.241 


0.0218753 
6.08769 
7.64991 
188.713 


0.0186789 
335.07 
28.0261 
218.586 


0.0253152 
357.142 
14.5679 
196.273 


0.0205366 
14.7651 
1.13126 
177.828 


TR (^) [GeV] 
\mrj 


1.77 X 10^1° 


2.63 X 10~^° 


2.50 X 10^1° 


4.02 X 10~^° 


7.3 X 10~" 


2.82 X 10^1° 



TABLE III. Best fit solutions for fermion masses and mixing obtained assuming the type-I seesaw 
dominance in the minimal SUSY 50(10) model. Pulls of various observables and predictions 
obtained at the minimum are shown for six different data sets. 



cases B, C, D with relatively large tan /3 but without inclusion of threshold correction 
give quite bad fit. There is a clear correlation between the overall fit and the presence 
or absence of the b-r unification in type-II models. Cases corresponding to the absence 
of the b-T unification cannot reproduce the atmospheric mixing angle and results in 
relatively poor fits. Inclusion of threshold corrections improves the fit but still ^ and 
the atmospheric mixing angle cannot be_ reproduced within la. The fit for tan/3 = 10 
obtained here with inputs from 



23, 



28 



is poor compared to the corresponding fit 



resented in [6|] which uses input from [21[. Compared to data in [21j, the result from 



|23[ | display larger deviation from the b-r unification and also errors in more recent 
input that we use for sin^ 6*33 are smaller. Both these features combine to give larger 
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pulls for the ratio ^ and sin^^js results in poor fit. 

• In contrast to the type-II case, the fits obtained in type-I case are uniformly better. 
Here one does not expect correlation between the atmospheric mixing angle and b-r 
unification. Thus the cases B, C, D with large tan /3 also give quite good fits. Even in 
these cases (except D) main contribution to comes from the pull in the atmospheric 
mixing angle. Threshold corrections are significant for large tan /3 and specific cases 
CI, C2 achieve b-T unification but the overall fit worsens compared to B, C, D. Unlike 
in the type-II case, the value obtained here for tan /3 = 10 is comparable to the 
corresponding value in p]. 

• We have fixed the overall scale of neutrino mass ri(r/j) in the case of type-II (type- 
I) seesaw by using the atmospheric scale as normalization. The resulting values are 
displayed in Table ( IITI IIIII) . ri(rR) arise from the vev of the components of 126// trans- 
forming as (3, 1, -2) ((1, 3, -2)) under the SU{2)l x SU{2)r x U{1)b-l- In particular, 
((1,3, — 2))y26^ sets the scale of the B — L breaking and is directly determined from 
the fits to fermion masses in the type-I scenario. From eq.(l9]), 

((1.3,-2))j26^ ^r^^i;s^cos/3 , (10) 

where Sm gives the mixing of the light Hd in the doublet part of 126// and v ~ 174 
GeV. r/j is roughly independent of the input data set and for the value r/j ^ 2.6 x 
10-^°m^/m2 GeV, eq.(IIOD gives 

((1,3,-2))^^ ^ 3.7 X lO^^s^ cos /3 GeV 

Thus the B — L breaking scale in the type-I seesaw can be close to the GUT scale for 
SmCos/3 ~ 0{1). It would however be significantly lower for large values of tan/? and 
would confiict with the constraints from the gauge coupling unification. The deter- 
mination of the B — L breaking scale in the type-II dominated scenario is dependent 
on the details of the superpotential. Earlier js], ^ analysis in the minimal model has 
shown that this scale cannot easily be lifted to the GUT scale and pauses a problem 
with the gauge coupling unification in the minimal scenario both for the type-I and 
type-II seesaw dominance jsl. Thus one does need to go beyond the minimal model 
and models with 120// are possible examples. 



B. Numerical Analysis: Extended model with 10 + 126 + 120 Higgs 

We now consider the non-minimal case obtained from eq.(|5]) by choosing all parameters in 
H, F, G as well as r, s, t^, t/j, r/j real. As before, the parameters r and r/j(r/,) determine 
the mt and overall scale of neutrino masses in type-I (type-II) seesaw dominated scenarios. 
Our choice of 14 observables is the same as in the previous subsection. But they are now 
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determined from the more general expression with non-zero G. H can be made diagonal 
without loss of generality. The mass matrices Mu, Md, Mi, Mb, Mr and Ml are expressed 
in terms of 16 real parameters (3 in if, 6 in F, 3 in G, s, t;, 1^-, and Id ) which determine 14 
observables Pi defined before. 





A 


B 


C 


D 


CI 


C2 


Observables 


Pulls obtained for best fit solution 


{mu/mc) 
{iTic/mt) 
{md/ms) 
{ms/rrib) 

(mb/nir) 


0.00196316 
0.000750815 
0.0547314 
0.0565403 
0.0114456 
-0.00279654 
-0.171035 


0.019005 
-0.114469 

0.618531 

0.473347 
-0.0155357 

-0.66999 
-0.301056 


-0.026015 
0.0964863 
0.0606721 

-1.30774 
0.0482971 

0.111235 
-0.381508 


-0.00109589 
0.296526 
-1.14305 
0.675173 
-0.00371258 
-0.0605957 
1.57926 


-0.00812155 
-0.0278823 
0.0271889 
0.0105556 
0.00167012 
0.00155096 
0.0514536 


0.0000225717 
-0.00413523 
0.00312586 
0.0361755 
0.000128709 
-0.00249006 
-0.048487 


sin 

sin 0^3 
sin0?3 
sin^ 6*^2 

sm (:^23 

5cKM [°] 


0.00833338 


-0.297416 


0.191515 


0.2129 


0.00050834 


-0.00412658 


-0.0106839 
-0.00295777 
-0.00466345 
0.0106277 
-0.0198083 
-0.00915099 


-0.0145213 
-0.058218 
-0.661544 
-0.194399 
1.08433 
0.168633 


-0.0420229 
0.301593 
0.381317 
0.333404 
-0.472589 
-0.520071 


0.0809603 
0.341191 

-0.632744 
0.399294 

-0.885401 
0.246618 


-0.00715584 
0.0120366 
-0.137308 
0.00217496 
0.0314484 
-0.0314877 


0.0000538731 
-0.000633901 
0.00650479 
-0.0043514 
0.00752103 
-0.0382519 


• 

A-rrwn 


0.0364 


2.9315 


2.7639 


5.921 


0.0254 


0.0038 


Observables 


Corresponding Predictions at GUT scale 


sin2 0i3 
(^MNsH 

«2n 

r^^jn,- [GeV] 


0.0215726 
34.3864 
6.26083 
161.011 

1.27 X 10"^ 


0.0312498 
5.21955 
76.0772 
253.288 
9.57 X 10"^° 


0.03568 

89.5 
289.921 
76.0613 
6.82 X 10"^° 


0.0214329 
315.898 
80.1968 
283.63 

1.56 X 10-^ 


0.0289663 
355.507 
60.3609 
220.306 

2.36 X 10^9 


0.0069694 
75.6953 
240.526 
34.4702 
3.68 X 10-9 



TABLE IV. Best fit solutions for fermion masses and mixing obtained assuming the type-II seesaw 
dominance in the non-minimal SUSY 50 (10) model. Pulls of various observables and predictions 
obtained at the minimum are shown for six different data sets. 

Results of numerical analysis carried out separately for the type-II and type-I dominated 
seesaw mechanisms are shown in Table (11 Vp and Table (jV|) respectively. The following re- 
marks are in order in connection with the results presented in these tables. As discovered in 
earlier numerical analysis jo], [isj, the introduction of the 120j^ leads to remarkable improve- 
ment in numerical fits in the type-II case. This mainly arises because the near maximality 
^23 is not directly connected to the the h-r unification. Thus the cases B, C, D which do not 
have the h-r unification also lead to very good fits in contrast to the minimal case. The fits 
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A 


B 


C 


D 


CI 


C2 


Observables 




Pulls obtained for best fit solut 


ion 




{mu/mc) 
{rric/mt) 
{md/ms) 
(ms/mh) 

{m^Jmr) 
{mb/nir) 


-0.0151499 
-0.000384519 
-0.0778857 
-0.052311 
0.00127584 
-0.0553488 
-0.0103881 


0.0262493 
0.000812518 
-0.0653974 

0.0706689 
0.00152407 
-0.0188764 

0.0214596 


-0.0019449 
0.000258262 
-0.0053692 

0.0726379 
0.00164957 
-0.0212797 

0.0260868 


-0.00461056 
0.00461629 
0.0272334 
0.0830354 
0.0268034 
-0.0282999 
0.0498589 


0.00542513 
-0.00304033 

0.00701785 

0.0120634 
-0.00722174 
-0.00866254 

0.00493891 


-0.0000775584 
0.0049013 
0.00147573 
0.0296307 
-0.000349497 
-0.00875371 
0.00967378 


( ^^lol \ 

sin 

sin Oin 
sin Of^ 
sin2 0^2 
sin^ 9^23 

5cKM [°] 


0.0324886 


0.00926157 


0.00312614 


-0.00630473 


-0.00302065 


0.000653399 


0.0159112 
0.0375281 
0.0309917 

0.00539037 
0.0332756 

-0.0585649 


-0.0140628 
-0.00674466 

0.0571306 
-0.0176765 

0.0143127 
-0.00882152 


-0.000195379 
-0.00216987 

0.175888 
0.00577816 
0.0125096 
-0.0406312 


0.00791696 
0.00501282 
0.0213394 
-0.013618 
0.0200216 
-0.0292954 


-0.0171517 
0.0100126 
-0.131639 
0.0092152 
0.00356131 
-0.0291351 


-0.000184021 
0.00590551 
-0.00184989 
0.000404734 
0.00026684 
-0.0310722 




0.0204 


0.0150 


0.0392 


0.0137 


0.0191 


0.0011 


Observables 


Corresponding Predictions at GUT scale 


sin2 

OMNSi J 

a2[°] 


0.0122064 
87.6747 
5.82048 
339.846 


0.0168745 
22.8731 
167.229 
331.88 


0.0146633 
330.351 
192.077 
34.5585 


0.0359278 
282.035 
286.062 
336.358 


0.0246489 
272.186 
352.828 
17.6585 


0.030277 
84.0238 
329.804 
325.182 


TR (^) [GeV] 
\mrj 


6.56 X 10^1^ 


1.22 X 10~i2 


1.34 X 10~^2 


3.03 X 10~i5 


5.0 X 10"^'' 


1.40 X 10~^3 



TABLE V. Best fit solutions for fermion masses and mixing obtained assuming the type-I seesaw 
dominance in the non-minimal SUSY SO{10) model. Pulls of various observables and predictions 
obtained at the minimum are shown for six different data sets. 

in cases (A, CI, C2) which have b-r unification are even better and all the observables are 
fitted almost exactly in these cases. These include the low tan /3 inputs and cases with large 
tan P and threshold corrections. As the results of Table (jVT) show, the fits obtained assuming 
the type-I seesaw dominance are uniformly better compared to the corresponding type-II 
results and show significantly improvement over the minimal model with type-I dominance, 
TablednH). 

One important difference compared to the minimal case is the overall B — L scale deter- 
mined from the neutrino masses. Unlike the minimal case, the values of in Table (jVT) are 
strongly dependent on the input data set and in some cases are quite large although each 
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data set appear to give very good fit to fermion masses. For example, one obtains in case 
(A) from eq.din]) and TabledV]), 

((1,3,-2))^^ ^ 1.5 X lO^Vcos/? GeV 

Thus reproducing neutrino masses in this case would require fine tuning Sm ~ 10~^ if the 
B — L breaking scale is to be close to Mqut- In contrast, in case C2 with tan/3 = 38, 
Table (jV]) gives 

((1,3,-2))^^ ^ 1.8 X IQ^hrn GeV 
which is close to the GUT scale. 



IV. FERMION MASSES IN NON-SUPERSYMMETRIC MODELS: NUMERICAL 
ANALYSIS 

We now turn to the discussion of various non-supersymmetric models. We shall consider 
three different cases. 



1. The minimal scenario containing the Higgs representations 45^ + 10/^ + 126//. 



2. Alternative model with 45 j/ + 120// + 126// proposed and analyzed in case of two 



generations in 26 



3. The non-minimal scenario with 45// + 10// + 126// + 120// with Hermitian structure. 

The case (2) is found to be unable to fit all the fermion masses and mixing angles. The 
minimal case works quite well in this regards and there is no real motivation to go to the 
non-minimal case as far as the fermion masses are concerned. We have included this for 
completeness and find that this case works even better than the minimal case. 



The 5*0(10) breaking [l7l] and the gauge coupling unification [16[ with intermediate scale 



has been reanalyzed recently following earlier works [18|, Il9| . The minimal case is argued 
to be adequate in achieving both the gauge coupling unification and breaking of 50(10) 
to the SM through an intermediate scale. The exact value of the required intermediate 



scale depend on the chain of the 5*0(10) breaking and various cases are given in [16[. The 
45// field contains components transforming as (15,1,1) and (1,1,3) under the Pati-Salam 
group 5f/(4) X SU{2)l x SU{2)r. They allow 50(10) braking to SU{?,)c x SU{2)l x 
SU{2)r X U{1)b-l, SU{A) X SU{2)l x U{1)r or to SU{b) x U{1) groups. At the tree 
level, only SU{b) intermediate stage is shown to lead to consistent spectra without tachyons 
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but this chain does not preserve the gauge coupling unification. As discussed in [Jj 
turning on 1-loop corrections also allows the other two breaking chains which preserve gauge 
coupling unification. The final breaking to SM can be achieved by the (1,1,3) component 
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of 126h which also leads to neutrino masses. The IOh and 126^ contain respectively bi- 
doublets (1,2,2) and (15,2,2). They need to mix in order to finally generate the standard 
model doublet(s) simultaneously containing the 10^;^ and 126// components. This can be 
achieved by fine tuning. For example, the mixing between bi-doublets is achieved through 
the following term 

y ~ XijXkl^ijklm4>m (11) 

which couples 45// (x) to 126// (S) and 10// (0). This mixes two bi-doublets when com- 
ponent of 45// transforming as singlet under the S't/(3)c x SU{2)l x SU{2)r x U{1)b-l 
acquires a vev. Then through fine tuning one can keep one of the two bi-doublets in 10// and 
126// at the intermediate scale. Subsequent breaking to SM is achieved through the (1, 1, 3) 
component of 126//. Eq. ffTTl) provides this way the required mixing between doublets in 
T26// and 10//. 



A. Numerical Analysis: Model with only 10 + 126 Higgs fields. 



A non-supersymmetric 5*0(10) model with 10 and 126 Higgs fields together with U{1)pq 
symmetry has the same Yukawa interactions as the minimal SUSY 5*0(10), eq.([5]) with 
G = 0. Minimization is performed based on the input values of the charged fermion masses 
obtained by running quark and lepton masses up to the GUT scale with m//=140 GeV 
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. We use the updated low energy values of quark mixing angles, CP phase and neutrino 
parameters since the effect of RG is known to be negligible for hierarchical neutrino spectrum. 
We reproduce all the input values in Table ( TVT|l for convenience. As before, we take and 
Ml as independent and express the remaining matrices in terms of them and r, s as in 
eq.([9]). Since the masses of the charged leptons are known precisely, we go to the basis with 
a diagonal Mi and use them as fixed input. Thus we have 15 real parameters (12 in M^, 
complex s and real r) which determine remaining 13 observables shown in Table (IVip . The 
function is defined in terms of these parameters. 



Results of numerical analysis carried out separately for type-I and type-II dominated 
seesaw mechanisms are shown in Table fIVIip . Parameters obtained for the best fit solutions 
are shown in Appendix A. It is evident that the type-II mechanism fails completely in 
reproducing the spectrum. Once again this is linked to the complete absence of the h-r 
unification in non-supersymmetric theories. Neither the atmospheric mixing nor the h quark 
mass can be reproduced correctly in this fit. In contrast, the type-I seesaw works quite well. 
In fact, the quality of fit in this case is much better than the minimal supersymmetric model 
with type-I seesaw. Table (lllip . 

As before the gets determined from the atmospheric neutrino mass scale. Assuming, 
that only one standard model survives at the electroweak scale one has, 

((1,3,-2))^^ ^r^hmV 
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GUT scale values with propagated uncertainty 


md(MeV) 


— U.4e5 


Ami, (eV2) 

sot V / 


(7.59 ±0.20) X 10-5 


mJMeV) 


22+^ 

— D 


Am2,^(eV2) 

(XVXTl \ / 


(2.51 ±0.12) X 10-3 


mb(GeV) 


1.00 ±0.04 


sin 


0.2246 ±0.0011 


m„(MeV) 


0.48+0-2? 

— U.l ( 


sin 


0.0420 ± 0.0013 


mc(GeV) 


n 9qc:+0.035 


sin 


0.0035 ± 0.0003 


mt(GeV) 


74.oi^:? 


sin2 


0.3208 ± 0.0164 


me(MeV) 


0.469652046 ± 0.000000041 


sin2^^3 


0.4t)29_o 0484 


m/,(MeV) 


99.1466226 ±0.0000089 


sin^ ^^3 


< 0.049(3(7) 


m^(GeV) 


1.68558 ± 0.00019 


^CKM 


69.63° ± 3.3° 



TABLE VI. Input values for quark and leptonic masses and mixing angles in the non- 
supersymmetric standard model extrapolated at Mqut = 2 x 10^^ GeV. 



tr in Table (IVIip gives 

((1, 3, -2))^^ ^ 3 X lO^^s^GeV . (12) 

Unlike the sup ersymmetric model, one would like to have this scale at an intermediate value 
10^^ GeV 16| in order to achieve the gauge coupling unification. This will require substan- 



tial fine tuning. The exact value of the required intermediate scale for the gauge coupling 



unification would depend on threshold effects not included in the analysis in 
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This would 



need a detailed study of the scalar potential minimization and the scalar sector of the theory. 



The leptonic parameters ^13 and three CP violating phases ai^2 and 5mns get fixed at 
the minimum and are shown in Table (lVIII) . The firm predictions on these observables in the 
scheme can be obtained by checking the variation of with the values of various observables. 
Following, , we pin down a specific value pq of an observable P by adding a term 

2 _ [ P-p, \' 

to and then minimizing 

-2 2 I 2 

X =X +Xp ■ 

If P happens to be one of the observables used in defining x^i then its contribution is removed 
from there. Artificially introduced small error fixes the value for P at the minimum of 
the x^. We then look at the variation of 

xl,in = i.Sd -Xl)\min (13) 

with pq. The results of such analysis carried out for the observables sin^ 6'23 and sin^ 9^^, 
are displayed in FigH] and Fig|2] respectively, sin^ ^23 can assume value in large range and 
the 90% confidence level bound corresponding to A^^ = 4.61 covers its entire 3(T range 
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Type-I 


Type-II 


Observables 


Fitted value 


pull 


Fitted value 


pull 


md 


0.000810163 


-0.687161 


0.00101285 


-0.264898 


Tils 


0.0208099 


-0.198354 


0.0225915 


0.0844982 


TUf, 


0.999667 


-0.00831657 


1.08201 


2.05031 


rriu 


00049.^02S 


0751 1 33 


000507336 


1 3668 




0.237348 


0.0670883 


0.237096 


0.0598882 


mt 


73.9427 


— (J. 1)154941 


74.3006 


U.(J7ol44 


TTT-e 


0.000469652 


— 


0.000469652 


— 




0.0991466 


— 


0.0991466 


— 


rrir 


1.68558 


— 


1.68558 


— 


[ ^^sol \ 
\'-^"''atm / 


0.030526 


0.127968 


0.0297114 


-0.235285 


sin Otc. 


0.224651 


0.0464044 


0.224499 


-0.0916848 


sin Bin 


0.0420499 


0.0392946 


0.0421308 


0.103004 


sin 0?q 


0.00349369 


-0.0974312 


0.00353053 


0.0389979 


sin^ e\o 


0.323245 


0.148134 


0.3108 


-0.610792 


sin^ 6*00 


0.435096 


-0.369178 


0.113306 


-7.02461 


sin^ OL 

io 


0.0244287 




0.0176863 




<5cKM [°] 


69.5262 


-0.0314447 


69.2051 


-0.128759 


5mns[°] 


318.465 




14.5386 






21.5053 




345.645 






215.128 




141.905 






5.62 X 10-14 




2.09 X 10-1° 








0.710777 




54.1197 



TABLE VII. Best fit solutions for fermion masses and mixing obtained assuming the type-I and 
type-II seesaw dominance in the minimal non-SUSY 50(10) model. Various observables and their 
pulls at the minimum are shown. All the masses shown are in GeV units. The bold faced quantities 
are predictions of the respective solutions. 

0.33 — 0.64. In contrast, a clear prediction emerges for the angle sin^ which preferentially 
lies in the range 0.015 — 0.03. 



B. Numerical Analysis: Model v^ith only 120 + 126 Higgs fields. 

We now consider an alternative model obtained by replacing 10^ with 120h in the min- 
imal model discussed before. This model is argued to be quite attractive and predictive 



when restricted to the second and the third generations 26|. It is thus interesting to see if 



the model works in more realistic case with three generations which require explanation of 
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0.35 0.40 0.45 0.50 0.55 0.60 0.65 



FIG. 1. Variation of Xmm with sin^ ^23 the minimal non-susy 50(10) model with Type-I seesaw. 



30 




FIG. 2. Variation of Xmin with sin^ in the minimal non-susy 50(10) model with Type-I seesaw. 



several new parameters. 



The fermion mass relations in this model are given by eq.([5]) with H = 0. F can be made 
real diagonal without loss of generality. Mi is not diagonal and the charged lepton masses 
are included in the function ([8]) unlike the previous case of the minimal model where 
they were set as input. Since the errors in the charged lepton masses are extremely small, 
the numerical optimization algorithm we use is unable to converge to the solution in finite 
time. Thus we set 10% error in charged lepton masses and minimize the with respect to 
16 (3 in F, 6 in G, real s and complex ti,tu,t£,) real parameters. 

Results of numerical analysis carried out separately for type-I and type-II dominated 
seesaw scenarios are shown in Table OVIIip . The detailed fits are quite different in two cases 
showing that a simple proportionality of the type-II and type-I contribution observed in the 
two generation study |26| does not hold in general. The model fails badly in reproducing 
the fermion mass spectrum in either case. Analytic study of the two generations lead in 
the model to a relation m.^ ~ Sm^. This is born out in the detailed numerical study with 
three generations as well. But this relation becomes one of the causes of the failure of the 



model as is clearly seen in the Table (IVIIip . Likewise, the numerical fits lead to nearly 
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Type-I 


Type-II 


Observables 


Fitted value 


pull 


Fitted value 


pull 


rrid 


0.000186192 


-1.9871 


0.000223284 


-1.90982 


nis 


0.00267758 


-3.2204 


0.00296063 


-3.17323 


nib 


0.844022 


-3.89946 


0.836471 


-4.08822 


rriu 


000480Q6 


004801 31 


00048341 2 


01 70595 




0.23454 


-0.0135291 


0.237869 


0.0819818 




'7/1 n r: o 

74.05o 


0.0132566 


73.891 


—0.0294532 


rrie 


0.000467656 


-0.0424983 


0.000475465 


0.123771 




0.0964545 


-0.271534 


0.101839 


0.271587 


rrir 


2.61149 


5.49314 


2.60147 


5.43367 


[ ^^sol \ 


0.0303749 


0.0605819 


8.59 X 10^^ 


-13.4841 


sin 9% 


0.224581 


-0.0172464 


0.224591 


-0.00790894 


sin 


0.0419722 


-0.0218756 


0.0420623 


0.0490417 


sin of. 


0.00354561 


0.0948516 


0.00353062 


0.0393252 


sin^ 


0.321216 


0.0243762 


0.320612 


-0.0124452 


sin2 


0.450311 


-0.0544896 


0.0375094 


-8.59228 


^CKM 


69.5526 


-0.0234639 


69.5794 


-0.0153481 






59.7934 




315.705 



TABLE VIII. Best fit solutions for fermion masses and mixing obtained assuming the type-I and 
type-II seesaw dominance in non-SUSY 5'O(10) model with 120+ 126 Higgs. All the masses shown 
are in GeV units. Various observables and their pulls at the minimum are shown. 



vanishing solar scale at the minimum in the type-II case. This becomes an added cause of 
very poor fits. It appears from the results that the renormalizable model with 45 + 120 + 126 
Higgs fields is not a good candidate to obtain even fermion mass spectrum in spite of its 
attractiveness at the two generation level 
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C. Numerical Analysis: Model with 10 + 126 + 120 Higgs fields. 

As before, we consider the case with Hermitian (Dirac) mass matrices. The mass relations 
are same as in eq.(l5]) with all parameters real. We have chosen the basis with a diagonal 
H. Ml is not diagonal in this basis and we parameterize it as Mi = UiDiUj with Ui being a 
general unitary matrix expressed in terms of three angles and six phases and Di is a diagonal 
matrix for the charged lepton masses. One can rewrite Mi in eq.(l5]) as 

3F - itiG = H- UiDiU} 
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Since F and G are real, the real and imaginary parts of the RHS separately determine F 
and tiG in terms of the charged lepton masses and parameters of H and Ui which are put 
back in eq.©. The remaining fermion mass matrices can be expressed in terms of 17 (3 in 
if, 9 in Ui, real r, s,ti,tu,tD) real parameters in the case of type-I seesaw dominance which 
determine 16 observables Pj shown in Table (IVip . One parameter to becomes irrelevant for 
the type-II seesaw case. We do the numerical analysis for this case and results are shown in 
Table (HXD. 





Type-I 


Type-II 


Observables 


Fitted value 


pull 


Fitted value 


pull 




0.00113968 


-0.000676838 


0.00108711 


-0.110189 


rus 


0.0219909 


-0.00150966 


0.0142689 


-1.28852 


rrib 


1. 


0.0000376219 


1.19665 


4.9162 


ruu 


0.000480133 


0.000666686 


0.000486627 


0.0331338 


rric 


0.235007 


0.000211758 


0.240819 


0.166268 


mt 


73.9997 


-0.0000888053 


77.4295 


0.857367 


rue 


0.000469652 





0.000469652 







0.0991466 





0.0991466 


0.220249 


rrir 


1.68558 





1.68558 


0.000124065 




0302402 


000545016 


0260106 


— 1 88556 


sm 


0.224601 


0.00105776 


0.224567 


-0.0304356 


sin 6*^3 


0.0420001 


0.0000431604 


0.0431393 


0.897068 


sin 


0.00351992 


-0.000308192 


0.00338234 


-0.509862 


sin2 


0.320821 


0.000292661 


0.278093 


-2.6052 


sin2 0^3 


0.453034 


0.000947066 


0.343286 


-2.26804 


sin^ 


0.0306736 




0.00538748 




(^CKM [°] 


69.6278 


-0.000660788 


72.7155 


0.935014 




355.719 




46.8148 






60.079 




60.6202 




«2n 


214.691 




250.978 






1.56 X 10-1^ 




3.43 X 10"^° 








- 10"^ 




44.0801 



TABLE IX. Best fit solutions for fermion masses and mixing obtained assuming the type-I and 
type-II seesaw dominance in the non-supersymmetric 50(10) model with 10 -|- 126 -|- 120 Higgs. 
Various observables and their pulls at the minimum are shown. All the masses shown are in GeV 
units. The bold faced quantities are predictions of the respective solutions. 

Parameters obtained for the best fit solutions in type-I case are shown in Appendix B. 
Unlike the supersymmetric case, the presence of 120// does not help in improving the fits 
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in the type-II seesaw dominated case. But the fits obtained for the type-I scenario are 
considerably better compared to the corresponding supersymmetric as well as the minimal 
non-supersymmetric case. Pulls in all observables are practically zero in this case. 



The predictions of the model for the observables sin^ 6*23 and sin^ are displayed in FigJ3] 
and FigJH respectively. Once again a clear predictionsin^ 6'i3 > 0.015 emerges in this case. 
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FIG. 3. Variation of Xmm with sin^ the extended model with Type-I seesaw. 
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FIG. 4. Variation of Xmm with sin^ in the extended model with Type-I seesaw. 



V. NUMERICAL ANALYSIS AND UNDERLYING FLAVOUR STRUCTURE 

Numerical analysis presented in the previous section has demonstrated viability of vari- 
ous S'O(IO) models in explaining the fermion masses and mixing. In the process, it has also 
provided us with specific structure of fermion mass matrices which can be used to obtain 
some insight into the underlying flavour structure. We discuss one specific case namely, the 
minimal non-supersymmetric model from this point of view. 

At the 5*0(10) level, Yukawa couplings H, F, G determine the flavour structure of various 
mass matrices. Thus any underlying flavour symmetry if it exists should get reflected in 
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the structure of these matrices. Specific structures for the Yukawa couphng matrices have 
been used to predict relations between the (hierarchical) quark masses and (small) quark 
mixing, see for example 
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30 



31 



In a large class of such models, the observed masses and 
mixing patterns among quarks are reproduced when elements of the quark mass matrices 
are expressed as powers of one or two expansion parameters. Following this, we try to look 
for a similar parameterization for the underlying matrices F, H in case of the minimal non- 
supersymmetric model. We choose the Cabibbo angle A = 0.2246 as a convenient parameter. 
Elements of F and H in this case are then found to have the following hierarchical structure 
in the basis with a diagonal Mf 



( 0.513ei-659'A4 



0.361e-i-257^A3 0.685e0-843'A2 
0.361e-i-257»A3 0.119eii43'A2 q_^^q^-2.12Mx 



-0.226i \4 



F = 0.278e 



2.561i 



GcV 



\ 0.685e"**'*-^*A^ 0.490e 
/ 0.802e 



0.470e 
\ 0.892e 



-2.123i 



A 



1 



A* 0.470e^-y'J'A'^ 0.892e 



-1.283j 



2.90i 



A3 2.359e0-5i5»A2 



-1.283i \2 



y 0.639e 



2.034i 



A 



0.639e2"34«A 
1 



(14) 



33 element turns out to be largest both for F and H and we have normalized other elements 
by its value in writing the above structure. Most coefficients in powers of A are roughly 
(9(1) except for the 22 elements. 



The above structure determined numerically here is suggestive of an underlying f/(l) 
symmetry used 3^ in the Froggatt Nielsen (FG) approach. Indeed a simple U{1) can 
explain the occurrence of various powers of A in eq. ffMj) . Consider a U{1) symmetry with 
the U{1) charges 2, 1, assigned respectively to three generations of 16i?-plet. Both 10// and 
126// are assumed neutral under this symmetry. In this case, the 33 elements of F, H arise 
from the renormalizable couplings lQ^F^^^F<t>H {4> = 10, 126). The 23 and 32 elements 
follow from the couplings 162f^Q3F<Ph^- Likewise, the two, three and four powers of rj 
respectively generate 0(A^, A^, A^) terms in eq. ([HD where A = ^, M being some underlying 
scale above the U{1) breaking scale (77) and rj is assumed to carry the U{1) charge —1. The 
quark mass matrices resulting from the above F, H also follow this simple pattern as in 
eq.dUD: 



= 0.9708e' 



0.6809i 



GeV 



/ 0.800ei-48i»A4 0.539e-i-503«A3 1.02Ae°-^^^'X^ \ 

0.539e-i-503*A3 0.699e2-204^A2 0.733e-2-37o»;^ 

V 1.024e°-597u2 0.7336-2-370.^ 1 

/ 0.6096 



Mu = 72.639e 



0.523i 



GeV 



1.585i\4 



-1.359i \3 



A-^ 0.796e 



0.741i 



0.4196-1-359^3 0.281ei-98i'A2 0.5706-2-225^^ 
\ 0.796eO-^4H^2 o.5706-2-225«A 1 
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This structure is already proposed and studied in [3l| as a possible explanation of quark 
and neutrino mixing and masses. Here it follows from a detailed analysis of this specific 
50(10) model. As shown in 3l|], such a form can reproduce the observed mixing and mass 
patterns for quarks. The expansion parameter chosen in 3l| is somewhat larger, A 



0.26. 



The Dirac neutrino mass matrix on the other hand is given by 



Md = 86.240e' 



0.210i 



GeV 



0_201e-0-959.A3 o.567e-0-222'A2 0.27ie-^-^^^'\ 
^ 0.382ei "°*A2 0.273e-i *25iA 1 y 



(15) 



The coefficients in front of various elements are anomalously small and thus M/j does not 
really share the same symmetry as the underlying Yukawa matrices. The Md and Mr ~ F 
conspire to produce a neutrino mass matrix which has an interesting form 



/ 1.339e2-543*;^3 o.878126e-0-662«A2 1.753ei-529'A \ 



= 0.087e 



-0.898'i 



rnr^ GcV 



„2.543i;^3 

0.878e-0-662*A2 

,1.529i^ 



\ 1.753e^ 



-0m2^y2 1.753el.529. 

0.800e~2.646« l.062e-i-4^^* 
1.062e-i-458* 1 



(16) 



Since we are working in a basis with a diagonal M;, the above matrix determines physical 
neutrino mixing and allows us to understand the leptonic mixing structure analytically. 
Firstly, the 23 block has all elements of 0{1) which results in the large atmospheric angle 
and hierarchy in neutrino masses. Secondly, the 11 and 12 elements are zero to leading order 
in A. In the approximation of neglecting higher powers of A, the has two-zero texture 



(classified as Al in 
mixing angle 33 



33 



The presence of the zeros leads to a firm prediction of the third 



sin^ e\ 



13 



AmL; \ sin^ 6*5^2 cos^ 6 



'12 



COS 26{2 tan^ ^23 



(17) 



This analytic relation is in very good agreement with the numerical values. Evaluation of 
the RHS using the best fit values of parameters in Table fl VIII ) leads to sin^ ^ 0.0245 in 
agreement with the numerical prediction. Even away from the minimum x^, one would get 
sin^^^g around 0.02 as long as two zero structure and hence eq. (fT711 holds approximately. 
This is born out quite well in Figure |2j 



The simple U{1) symmetry used to explain the structure of F, H may appear to have 
two shortcomings. Firstly, the specific structures are found in a basis with a diagonal Mi 
Secondly, the coefficients of powers of A in F, if are not strictly 0(1), notably in the 22 
elements. In general, the definition of symmetry and resulting texture of Yukawa matrices 
are basis dependent. Basis with a diagonal Mi are very special basis and it would be more 
desirable to find a basis in which Mi also has a structure similar to the F, if, M^, M„. One 
can indeed find a class of unitary rotations which bring the diagonal Mi to the form as in 
eq. (fT^ and at the same time retain the forms of F, H albeit with a different set of coefficients. 
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The U{1) symmetry leads to the following general form of the Yukawa matrices: 



h3 



F,H 
hi 



12 



V 



F,H , 



13 
F,H 



a. 



23 



A 



(18) 



''23 



F, H as given above can be diagonalized with high accuracy by rotation Rf,h consisting of 
three successive rotations in 2 — 3, 1 — 3 and 1 — 2 plane with the mixing angles 31 



sin 6. 
sin 6 

tan 261 



F,H 
23 
F,H 
13 

F,H 
12 



F,H ^ 
^23 ^ 



F,H ^2 
'^13 ^ ' 

2A 



F,H 
^12 



F,H F,H 
^23 '^13 



F,H 
^22 



— a: 



^.^^^2 

23 / 



0(A2 



(19) 



The eigenvalues of F,H are (9(1,A^,A^). The eigenvalues of Mi are roughly of similar 
order-though the coefficient for the electron mass is somewhat small. Thus Mi can be put 
to the form as in f[T^ by rotating the diagonal Mi with a rotation matrix Vi with angles as 
in eq. f|T9|) but with a different set of coefficients a[j. It is easy to see that when F,H are 
expressed in new basis their forms do not change to leading order in A but now coefficients in 
front of powers of A are different say, CLif^- They depend on a^'^ and a\j. Thus symmetry in 
question may manifest itself in more general basis than the specific diagonal basis provided 
a-J'^ are also C(l). 



Let us consider a simple example. Rotate F, H and diagonal Mi with a common rotation 
Vi defined as 

/l \ 



Vi 



1 

- 



- 1/2A2 

-Xe~'l^ 1 



Xe'P 

- l/2X^ I 

P can be chosen such that the coefficient of various powers of A in elements of F' = Vi^FVi 
and H' = ViHVi are near to 1. The best fit value of (5 turns out to be (3 
this one gets 



1.055 and for 





( 0.802A'' 


0.772A3 


0.892A2 


\ 


F'\ = 0.278 GcV 


0.772A3 


1.12A2 


1.638 A 






^ 0.892A2 


1.638A 


1 


/ 




( 0.513A4 


0.593A3 


0.685A2 


\ 


H'\ = 1.088 GcV 


0.593A3 


0.939A2 


0.876A 






^ 0.685A2 


0.876A 


1 





(20) 



Unlike in eq. (ll4l) . all the cofficients of various elementsin the above equation are now 0{1). 
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Ml is non-diagonal in this basis and is given by 

/ O.IOQA^ 



\M,\ = 1.685GeV 



v 









1.06A2 1.006A 
1.006A 1 



(21) 



We note that the Yukawa couphng matrices in cases other than the minimal also display 
hierarchical structure, see results in Appendix (B). Thus these cases can also be understood 
in terms of some simple pattern as in the minimal case discussed here explicitly. 



VI. SUMMARY 



S'O(IO) models have been used to obtain a unified description of fermion masses and 
mixing angles. We have undertaken in this paper an exhaustive analysis of many different 
SO (10) models. Using several different data sets as input, we have numerically determined 
viability of these models in reproducing the fermion spectrum. In case of the supersym- 
metric models, we used data corresponding to different values of tan (3 and with or without 
appreciable finite threshold correction. Comparison of different set clearly brings out an 
important feature. In the minimal model with type-II seesaw dominance, the h-r unification 
appears to be a key ingredient. The cases without such unification cannot explain the entire 
fermion spectrum. In particular, the case of very low tan /3 showing this unification works 
much better than the previously studied data set with tan/3 = 10. This connection is not 
required if neutrinos obtain their masses from the type-I seesaw mechanism. In this case 
one can obtain very good fits in the minimal model almost for every data set used, see 
Tabel(III). Moreover, the B — L breaking scale inferred from neutrino masses also lies closer 
to the GUT scale compared to the type-II seesaw mechanism. The situation becomes better 
when a 120-plet of Higgs field is added to the model. Here one can get excellent fits to 
fermion masses in both the type-I and type-II seesaw mechanisms. 



We also carried out a detailed analysis of the fermion masses in non-supersymmetric mod- 
els. The minimal non-supersymmetric model with 45// + 10// -|- 126// is quite economical and 
is argued recently [16|, ll7| to be a viable candidate for the gauge coupling unification. As 
shown here it also provides a very good description of fermion masses as well. Intermediate 
scale ~ 10^^ GeV is required in this model in order to obtain the unification of gauge coupling 



16[. The scale preferred from the fits to fermion masses presented here is somewhat larger. 
This scale can be reduced if the admixture of the light doublet in the doublet component 
of 126// is very small, see eq. (fT2|) . Viability of these as well as simultaneous analysis of 
the constraint from the gauge coupling unification will depend on the detailed analysis of 
the scalar sector of the theory. The Yukawa coupling matrices obtained numerically in this 
case display interesting structure which can be understood from a very simple symmetry 
imposed at a high scale. These features coupled with its economy makes the minimal non 
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supersymmetric model an attractive choice to unify basic gauge and Yukawa interactions. 
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VII. APPENDIX 

We list here the fermion mass matrices using the best fit values of the parameters given in 
Table IVIII (Table lIXp corresponding to the type-I seesaw mechanism in the case of minimal 
(non-minimal) non-supersymmetric 5*0(10) model. All the mass matrices are expressed in 
GeV units. 



A. Best fit parameter values: The minimal nonsusy SO{10) model, type-I seesaw 
mechanism (Table [Vll|) . 

Parameters obtained for best fit solution. 



r = 69.1739; s = 0.362941 - 0.0463175i 
/ 0.000469652 \ 

Ml = 0.0991466 

\ 1.68558 y 

/ -0.00110182 + 0.00164125i 0.0040374 - 0.00434507i 0.0145011 + 0.0480084i \ 



M, 



0.0040374 - 0.00434507i 
0.0145011 + 0.0480084i 



-0.0331074 + 0.00870484i -0.0187112 - 0.158707i 
-0.0187112 - 0.158707i 0.754282 + 0.611126i / 



(22) 



Results: 



/ -0.0575896 + 0.0967087i 0.231322 - 0.25593i 0.881792 + 2.78041i \ 

0.231322 - 0.25593i -0.826159 + 0.612181i -1.21531 - 9.21493i 
^ 0.881792 + 2. 78041i -1.21531 - 9.21493i 62.9262 + 36. 2872i y 

/ -0.0000981928 + 0.00131563i 0.0000421662 - 0.003853i 0.0276444 + 0.0202258i \ 

-0.0640659 + 0.0272358i -0.0653091 - 0.0653272i 
-0.0653091 - 0.0653272i 0.054234 - 0.0680089i j 



V 



0.0000421662 - 0.003853i 
0.0276444 + 0.0202258i 
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B. Best fit parameter values: The non-minimal nonsusy 50(10) model, type-I 
seesaw mechanism f Table IIX|) . 



Parameters obtained for best fit solution. 



H = 



F = 



G = 



-52.4173; s = 1.61949; U = 3.1751; <„ = 0.0413014; to = -11.7339. 
/ 0.00158452 ^ \ 

0.0407501 
\ -0.330398 y 

/ -0.00116221 -0.000145513 0.0130876 \ 
-0.000145513 -0.0224155 -0.00121344 
0.0130876 -0.00121344 -0.667509 / 

( -0.00670763 0.00612927 \ 

0.00670763 -0.0437162 

\ -0.00612927 0.0437162 



(23) 



Results: 



M„ = 



Ml 



/ 0.00042231 -0.000145513 - 0.00670763i 0.0130876 + 0.00612927i \ 

-0.000145513 + 0.00670763i 0.0183346 -0.00121344 - 0.0437162i 

0.0130876- 0.00612927i -0.00121344 + 0.0437162i -0.997907 y 

/ 0.0156028 0.0123525 + 0.0145214i -1.11099 - 0.0132693A 

0.0123525 - 0.0145214i -0.233181 0.103008 + 0.0946413i 

-1.11099 + 0.0132693i 0.103008 - 0.0946413i 73.9827 y 

/ 0.00507117 0.00043654- 0.0212974i -0.0392628 + 0.019461 A 

0.00043654 + 0.0212974i 0.107997 0.00364033 - 0.138803i 

-0.0392628 - 0.019461i 0.00364033 + 0.138803i 1.67213 y 

/ 0.243898 + 0.00702837i 0.0907917 - 0.0237474i -1.65214 - 0.184893i \ 
r^r^ 0.0907917- 0.0237474i 4.6052 - 0.0433382i -5.76376 + 0.720305i (24) 
-1.65214- 0.184893i -5.76376 + 0.720305i 5.33677 + 1.11414i y 
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